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Abstract 

Kernel functions related to quantum many-body systems of Calogero-Sutherland 
type are discussed, in particular for the elliptic case. The main result is an ellip- 
tic generalization of an identity due to Sen that is a source for many such kernel 
functions. Applications are given, including simple exact eigenfunctions and cor- 
responding eigenvalues of Chalykh-Feigin-Veselov-Sergeev-type deformations of the 
elliptic Calogero-Sutherland model for special parameter values. 
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1 Introduction 

Sen proved in [Ij that the quantum many body Hamiltonian 



^:=-E;^^+ E 1JkV{Xj-X,) (1) 

J=l "^J 1<J<K<M 

for V{r-) = l/[4 sin^(r/2)] and the coupling constants 



IJK ■= A(mj + mK){XmjmK - 1) (2) 
has the following exact groundstate 

$o(X):= J] e[Xj-XKf^-^^^ (3) 

\<J<K<N 

with 6'(r) = sin(r/2), for arbitrary particle number A/", coupling parameter A > 0, and 
particle masses mj/2 > 0. He was also able to compute the corresponding ground state 
energy So exactly pLj . Sen's result is a generalization of a well-known result for the Calogero- 
Sutherland (CS) [2113] model (corresponding to the special cases where the particles are non- 
distinguishable and nij = 1 for all J). As discussed in [1], Sen's identity {'H — Sq)^o(X.) = 



^ langniann@kth.se 
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holds true for arbitrary real parameters mj and A, and this (obvious) generalization is inter- 
esting since it has many special cases that provide means to compute exact eigenfunctions 
and corresponding eigenvalues of the CS model and, more generally, of Chalykh-Feigin- 
Veselov-Sergeev- (CFSV-) type deformations [SH7] of the CS model. The latter are given 
by differential operators 

j=l i l<j<k<N 

J=^ 1<J<K<N 
N N 

j=l K=l 

and reduce to the CS Hamiltonian for = (A^ and are non-negative integers such 

that + > 1, and we write x short for (xi, . . . ,xn) and similarly for x). The above 
mentioned special cases are identities of the form 

(if^^^(x, i) - if,,,M(y, y) - C) ^(x, i, y, y) = (5) 

for some function F and some constant C. We refer to such F as kernel function of the 
pair of differential operators (if^^(x, x), /f^^ jQ^(y, y)) in the following. 

In this paper we give a generalization of Sen's identity to the elliptic case, i.e. to the case 
where V{r) is (essentially) the Weierstrass elliptic function p(r). We also point out various 
special cases of this elliptic identity that provide means to compute exact eigenfunctions 
and eigenvalues of the CFSV-deformed elliptic CS (eCS) differential operator in (j4]), in 
generalization of results in |4j (trigonometric limit) and [HHTT] (eCS model). 

Our Corollary 12.21 generalizes and unifies identities given in [12], (6)-(12) (the latter 
correspond to the special cases N = M = and N = M = 0). Generalizations of the 
identities in [12] to the relativistic (Ruijsenaars) generalization of the CS-model [13j with 
a restriction on parameters as in 0161) were recently given by Komori, Noumi, and Shiraishi 
[H]; see also [15] and references therein for related results. Identities like in ([5]) were also 
used in other works to construct eigenfunctions CS-type operators, including [T ^ [T6 | [T7]. 

We find in this paper exact eigenfunctions of elliptic CS-type differential operators that 
can be represented by simple, explicit formulas (Corollaries 13.11 and 13.21) . The only other 
similar result in the literature we are aware of are BC^ variants of the eCS model [19] 
(also known as Inozemtsev model [18j) that were found by Gomez-UUate, Gonzalez-Lopez, 
and Rodriguez to be quasi-exactly solvable for certain special parameter values [20]; see 
also [HIES]. 

The results in the present paper and in [14J suggest that there should exist a generaliza- 
tion of Sen's identity to the relativistic case, and this identity might be easier to prove than 
its special cases found in [H]. This result would also be interesting since, as we expect, 
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(4) 



it should provide means to obtain results for CFSV-type deformations of the Ruijsenaars 
models. We also conjecture that there exists a Sen-like identity for Inozemtsev-type models. 

The rest of this paper is organized as follows. In Section |2] we present our results, and 
in Section |3] we discuss applications. Proofs can be found in Section HI Some results on 
elliptic functions that we need are collected in an appendix. 



2 Results 

We find it convenient to fix the elliptic function periods as 2uJi = 27t and 2uj2 = i/3 with 
(3 > 0, and we add a constant to p{r) so that the trigonometric limit /3 — )■ oo becomes 
simpleo 

^(^) ■= E A • l-R W91 = ^^^1^' + ^0 (6) 

4 sm [(r + ipmj /2\ 



with 



We also need 



_ ^1 _ ^ _ ^ 

Co :- — - ci - 2^ ^ . , 2.._ :- — . {7) 



^^2sinh^(/3m/2)' 12 



e{r) := sin(r/2) Y[{1 - 2g''" cos(r) + g^™), q := e"^/' (8) 

m=l 

proportional to the Jacobi theta function t9i(r/2, q). We use the short hand notation 



m 



:=^m}, n = l,2,3. (9) 



j=i 



Our main result is the following. 

Proposition 2.1. Let JV E N, and A and mj complex and non-zero for J = 1,2,... , A/". 

Then the differential operator in ([2])^([^ with V{r) in ^ and the function $o(X) in ^ 
with 6{r) in ^ obey the identity 

H + 2A|m|^ -foVo(X) = (10) 



'd/3 



with the constant 



£o = A2(^(|m^||m| - |m^|)co + (|m|2 - |m2|)|m|ci) (11) 

and Co,i in 

^Our conventions for special functions are as in [23]. Details can be found in [10^, Appendix A. 
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(Proof in Section UT]) 

In the following we are mainly interested in cases where A > and all rrij are real. 

Due to the /3-derivative term in f lTOj) it is not possible to interpret $o(X) as eigenfunc- 
tion of "H in general. However, regarding r := /3/(2|m|A) as (imaginary) time, one can 
interpret ( ITOj) as (Wick rotated) time evolution equation for quantum many-body systems 
with peculiar time dependent two-body potentials, and Proposition 12.11 gives an exphcit 
solution of this equation. It could be interested to explore this interpretation of (ITOj) and 
its consequences further. However, in the rest of this paper we only discuss a very different 
application: one can use (fTOj) as a source to obtain kernel functions for pairs of eCS-type 
differential operators as in (jl]) [1]. 

The key observation is this: One can choose the parameters and rrij such that "H 
in (II]) is (essentially) a difference of two operators if^jy(x, x) and iJj^^ ^^(y, y) Indeed, 
7ji<: = in the following four cases: (i) mj = 1 and rriK = —1, (ii) frij = 1 and rrix = 1/A, 
(iii) rrij = —1/A and tuk = —1, and (iv) rrij = —1/A and tjik = 1/A. One thus can divide 
the variables X in four groups x, x, y and y where the "mass parameter" mj is the same for 
all variables in each group. Choosing the mj as 1, —1/A, —1, and 1/A in these four groups 
the first two groups of variables decouple from the last two groups, and Proposition 12.11 
implies the following. 

Corollary 2.2. Let N, N, M, and M be non-negative integers and A complex and non-zero. 
Then the function 

FN,N,M,Mi^' ^' y' y) •=^o''^(^' x)^o^'^^(y, y) 

(n7=i nil oix, - m)) (nil utu oixj - y,)) (12) 



X 



(n7=i utU oix, - y,v) (nil nil oixj - m: 



1/A 



with 



v^^'^(x, i) := ^ (13) 



Ul<j<k<N ^(Xj - XkV^ (Ul<J<K<N ^{Xj - Xk 

uf=i nil ^ixj - xj) 

and the differential operators if^jy(x, x) and if^^ jQ^(y, y) in obey the identity 

d_ 

'dp 



(if^_^ (x, i) - H^ ,^, (y , y ) + 2 [(iV - M) A - iV + M] — - C^,M,^,Af ) F^,N,M,M (x, ^ , y , y ) = 

(14) 



with the constant 

i2 



Cn,m,n,m = [[N{N - 1) - M{M - 1)]X' - (N + M){N - M)\ + {N - M){N + M) 

- [N{N - 1) - M(M - 1)]/a)co + ((iV - M)[{N - Mf - N - M]\^ 

- [3(iV -Mf -N - M]{N - M)A + (iV - M)[3(iV ~ Mf - N - M] 

-{N -M) [{N -Mf -N - M]/)^ ci (15) 

and Co,i in 



(Proof in Section 

In cases where there is no /3-derivative terms, i.e. if the parameters are such that 

{N-M)\ = N-M (16) 

holds true, we obtain an identity as in ([S]), i.e. Fj^^ j^^jj^^ in f[T^ is a kernel function for 
the pair of differential operators (ifjY^(x, x), Hj^,^ j^j{y, y)) if flTB]) holds true. Note that the 
restriction in (fT6!) is not present in the trigonometric limit |4j. 

Remark 2.3. One can show that by redefining the elliptic functions 

V{r) p{r\7r,il3), e{r) ^ ^{r/2,q) (17) 

the identities in Proposition \2.1\ and Corollary \2/A hold true as they stand hut with the 
constants redefined as follows 

£q X(Af- l)|m|co 

^ ^ ^' ' ^ - - (18) 

Cj^^f,^M^M^iN + N + M + M-l)[\{N-M)-N + M]co 

(see Section \4 . 3\ for details). This shows that choosing standard elliptic functions makes 
the constants in our result significantly simpler. Our choice has the advantage that the 
trigonometric limit /3 — )■ oo zs obvious, whereas standard elliptic function requires a non- 
trivial multiplicative renormalization. Moreover, for the cases of interest to us the constants 
are simple anyway; see the remark below. 

Remark 2.4. We are mainly interested in the special cases of / fiOj) where |m| = 0. In such 
a case the constants in ( [77]) and simplify significantly as follows, 

£o = — A^|m^|co 

Cn,n,m,m = [-A'(iV -M) + {N- M)/X]co. ^^^^ 

Remark 2.5. At first sight it seems one could obtain a generalization of the identity in [I4\ ) 
from Proposition \2.1\ by choosing the "mass parameters" for the four groups of variables as 
m, —l/{mX), —m, l/{mX) with arbitrary real m ^ 0. However, multiplying the identity 
thus obtained by m and changing A to A/m^ one recovers [14\ )- We thus set m = 1 without 
loss of generality. 

Remark 2.6. Replacing the set of parameters {N, N, M, M, A) by {N, N, M, M, 1/A) leaves 
the identity in [I4\ ) invariant. This suggests that the duality transformation {N, N, A) — )■ 
(A^, A^, 1/A) should be an interesting symmetry of the differential operators in This 
symmetry is a generalization of the well-known duality of the Jack polynomials (see e.g. 
corresponding to the special cases /3 — )■ oo and N = 0; see /^/ for further details and results 
in the trigonometric limit. 

In applications it is convenient to use the following slight generalization of the result in 
Corollary E^l 
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Lemma 2.7. The identity in remains true as it stands if one replaces 

^ N,N,M,M\-^j -^^ y J y J ^ '^■^ N,N,M,MK-^^ -^^ y ^ y )^ 

Cn,m,n,m ^ Cj,^M,N,M + [N-M-{N- M)/Xy 
with arbitrary constants v eM. and c G C \ {0}, with |x| := Y2f=i^j 



(20) 



(Proof and physical interpretation in Section 

As will become clear in the next section, this result allows one to remove exponential 
factors, corresponding to trivial center-of-mass contributions, from eigenfunctions. 



3 Applications 

In this section we point out various interesting special cases of our results in the previous 
section. 

As already mentioned, it is not possible in general to interpret f IlOp as an eigenvalue 
equation. However, in case the parameters are such that |m| = this is possible. One 
such case of interest to us is obtained from Corollary 12.21 setting M = M = and denoting 

^N,Nfl,0 Eq-. 

Corollary 3.1. Let N and N he non-negative integers and 

X = N/N. (21) 

Then \&^'^(x, x) in /[T^) with 6{r) in ^ is an exact eigenfunction of the differential operator 
if^^(x, x) in 0) with V{r) in and the corresponding eigenvalue is 

Eo = {N- N^/N)cq (22) 

with Co in 

Another identity where one can construct simple eigenfunctions of the differential oper- 
ator H^ f;{x, Si) is ([11]) for M = 1, M = 0, and (A^ - 1)A - iV = 0, i.e. 

(i/^,^(x,x) + ^ - C^,^_,,o)cvl/^'^(x,5)P(x,5,y)e'''(''^l-^-I^IA) (23) 

with C^,^,i,o = [-A^(Ar - 1) - N/X]co = [A^ - 1 - Ny{N - l)]co and 

N N 

P(x, i, y) = (J[ e{x, - y)-') (n - y)) (24) 
j=i j=i 
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with constants w G M and c G C \ {0} to be determined; we used the generahzation of (fT^ 
pointed out in Lemma [2.71 Inserting 9(x) = (i/2)e~^^/^^(e'^) with 

oo 

e{z) := (1 - z) JJ (1 - q^"'z){l - g^™/^) (25) 

m=l 

analytical and non-zero for < \z\ < 1, we observe that the function in (!24l) is equal 
to P(z,z,0 := UjOizj/O-^UjOizj/0 times exp(i[A(|x| - y) - |i|]/2) with ^ := up 
to some finite and non-zero multiplicative constant; here and in the following we use the 
notation 

Zj:=e'''^, z:= {zi,...,zn) (26) 

and similarly for x (we inserted A^A — = A to simplify the exponential factor). We 
thus find that, for v = —A/2 and suitable c 7^ 0, c'P(x, x, y) exp (if [|x| — y — |x|/A] in 
(125]) is identical with 'P(z,z,^). This latter function can be analytically continued to the 
annulus 1 < |.^| < in the complex ^-plane and is there equal to its Taylor series 

J2nez^~^'^ri{'^,i) with functions P„(z,z) that can be computed by contour integrals; see 
f l29|) below. Inserting this and [d"^ /dy'^)^~"' = —n'^C,'"' and comparing equal powers of ^ we 
obtain the following. 

Corollary 3.2. Let N > 2 and N > 1 be integers and 

X = N/{N-1). (27) 

Then the differential operator if^^(x, x) in @) with V{r) in ^ has the following exact 
eigenfunctions labeled by integers n, 

^„(x,i) = <'^(x,i)P„(z,z) (28) 

with \E'^'^(x, x) in [Wjj and 6{r) in and 

^n(z,s) = j ^r(n^(e'^vo-')(n^(^'"/o) (29) 

with 6{z) in and the integration contour a circle \^\ = R of radius 1 < R < 
Moreover, the corresponding eigenvalue is 

E{n) =n^ + [N -I- N'^/{N - l)]co (30) 

with the constant cq in (0). 

There are many generalizations of the results above that are, however, more complicated 
in general. A particularly interesting special case of the identity in ( 1141) is for N = N and 
M = M, i.e. 

(^7v,7v(x,i) -^7v,7v(y.y))^7v,7v(x,i,y,y) = o (3i) 

with Fj^ := -F^v TV TV TV tjl^j) . The latter generalizes a result in [8j for the eCS model (spe- 
cial case N = 0) that can be used to construct a perturbative solution of the eCS model 
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to all orders [TOl[TT]. It should be straightforward to generalize this solution and construct 
eigenfunctions and corresponding eigenvalues of the CFSV-deformed eCS differential oper- 
ators. Similarly as in |3], one should also be able to construct different representations of 
these eigenfunctions and eigenvalues starting from any identity in f[T^ whenever (ITB]) holds 
true. The complexity of such a representation is determined by M + M: the smaller the 
latter the smaller the complexity The results in Corollaries 13.11 and 13.21 correspond to 
cases with the smallest possible complexities M + M = and 1, respectively. It would 
be interesting to study these solutions in more detail, but this is beyond the scope of the 
present paper. 



4 Proofs 



This section contains the proofs of the results stated in Section |2l 



4.1 Proposition 12.11 

The proof of Proposition 12.11 below is a straightforward computations using functional 
identities of elliptic functions collected in Appendix |Al 

With $0 in © we compute 

1 1 (9^ 

^■=^Y.—l^^o (32) 

where we suppress the common argument X of the functions W and $0; here and in the 
following. Straightforward computations give 

W = XmK(l>'iXj -Xk) + Y1 >^rnjmK<PiXj - Xk) J2 - ^l)) (33) 

J=l K^J K^J L^J 

with 0(r) in (jH]) and 0'(r) := d(f){r)/dr; we used We write W = Wi + W2 with 

Wi = ^ (^\{mj + mK)(p\Xj - Xk) + X^mjuiKimj + mK)(pixj - xx f^ (34) 

l<K<J<Af 

the sum of all two-body terms and 

W2 = Yl X'mjmKmL {<p{Xj - Xk)<P{Xj - X^) 

1<J<K<L<JV (^35) 

+ <j){XK - XMXk - Xj) + <P{Xl - Xj^Xl - Xk)) 

all three-body terms; all sums were made symmetric with respect to the summation indices 
using (|52|) . Inserting the identities in fHHj) and fj49|) in fl34|) we obtain 

m= {lJKV{Xj-XK)-X'mjmK{mj + mK)[2f{xj-XK) + Co]) (36) 

l<J<K<Af 
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(37) 



with in ([2]). Using ( !5T!) and (152|) for x = xj — x^, y = — (a^j — xl) and 2; = Xi^- — xl 
we find 

>V2 = - ^ \^mjmKmL[f{xj-XK) + f{xj-XL) + f{xK-XL)] 

1<J<K<L<M 

= -2 ^ ^ \^mjmKmLf{xj - Xk)- 

1<J<K<M Li^J.K 

Adding Wi and W2 and inserting ^YliLi^j k'^^^ ~ 1^1 ~ "^J ~ gives 

>V= ^ (^7jK^(Xj-A:x)-2A^|m|mjmi^/(xj-XA')-A^Comjmx(mj + mA-)). (38) 

i<j<x<Ar 

From ([3]) and (l50l) we conclude 

1 9 

^^"^"0 = - 5Z Amjm7^[/(xj - Xi^) - Ci] (39) 
and thus (132!) and (!38|) imply 

E 7..V(A-,-A-..) + 2A|m|i-|4„-£. (40) 

^ J=l l<J<_K'<Ar ^ 

with the constant 

£0= ^A^mjm;^(mj + m/^)co + 2|m|A^mjmi^Ci^ . (41) 

l<J<ft'<A/' 

This is equivalent to ffTOj) with in ([1]) and in f lTT]) - Q 



4.2 Corollary [2721 

An outline of how to obtain Corollary 12.21 from Proposition 12.11 is given in the main text. 
For completeness we provide the formal details here. 

We choose A/'=A^ + iV + M + M and set mj = 1 for 1 < J < A^, mj = -1/A for 1 < 
J-N <N,mj = -1 for 1 < J-N-N < M, and mj = 1/A for 1 < J-N-N-M < M. 
Denoting Xj as x j for 1 < J < A^, xj for 1 < J - A^ < iV, ?/j for 1 < J - A^ - A^ < M, and 
ijj for 1<J — A^ — A^ — M<Mwe find by straightforward computations that in ([1]) 
and ([2]) is equal to if^^(x, x) — Hj^j jy^{y,y) as defined in (jl]), $o(X) in ([3]) is proportional 
to -FAr,7v,M,A?(^' y' f) (12) 5 and So in is equal to C^^at^m.m (ffSD- Since obviously 
2A|m| = 2[A(A^ - M) - N + M] this implies the result. □ 



4.3 Remark [2731 

Equations ([1]), ([3]) and ffTOl) imply that, redefining the elliptic functions by /9-dependent 
constants 

V{r) V{r) + bo 

e{r) B^9{r), h := T^^logSi ^ ' 



changes the constant in (HT!) as follows, 

£q ^£o - ^ 'jj,Kbo - 2A|m| ^ XmjmKh 

J<K J<K (^43-) 

=So - (A'(|m2||m| - |m=^|) - A(Ar - l)|m|)6o - A'(|mp - |m2|)|m|6i 

(we inserted ([2]) and made a straightforward computation). The claim made in Remark 12.31 
corresponds to the special case 60 = co and 61 = ci (some further details are explained 
in [12], Remark 1.1). □ 



4.4 Lemma 12.71 

We note that the identity in (11 01) is invariant under 



$o(X) ^ c$o(X)e^^^^=i'"-'^^ 
£0 ^ £0 + |m|f ^ 



(44) 



for real v and non-zero complex c. Indeed, the function $o(X) in ([3]) is obviously invariant 
under translations Xj — )■ Xj + a, for all a G M, and this implies 

i:^$o(X) = (45) 
j=i 

and 

?^$o(X)e'^^^=i"^^^^ = e'''^^=^'^'^'(n + |m|t;2)$o(X). (46) 

Using the generalization of flTUl) obtained by the substitutions in (jH]) and restricting to the 
special case as in Section 14.21 we obtain the result in Lemma 12.71 □ 

The substitution in fj44l) has a natural physical interpretation as a change of the center- 
of-mass velocity of the particle system. 

Note added: The integrability of the CFSV type deformation of the eCS model in (jl]) for 
J\f = 1 and V{r) in (jH]) was proved in [25jFI 

Acknowledgments. I am grateful to Martin Hallnas for helpful comments. This work 
was supported by the Swedish Science Research Council (VR) and the Goran Gustafsson 
Foundation. 



A Elliptic functions 

For the convenience of the reader we collect here various properties of the elliptic functions 
V{r) in ([6]) and 6{r) in ([8]) that we use in Section [2l 

•^I thank an anonymous referee for pointing this out to me. 
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The function 
obeys the relations 
and 

with the function 



0(r):=|-log(0(r)) (47) 

|-0(r) = -V{r) (48) 

<p{rf = V{r) - 2/(r) - cq (49) 

/(r):=-^log(^(r)) + Ci (50) 
and the constants Co,i in ([7]). Moreover, 

0(x)0(y) + 0(x)0(^) + <P{y)(t>{z) = fix) + f{y) + /(;2) if x + y + ;2 = (51) 
with the same function / in ( 15U]) . We also use 

0(-r) = -0(r), V(-r) = V(r), /(-r) = /(r). (52) 
All these identities are classic; see [TU] and [12] for references and elementary proofs. 
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